Moment of Inertia and Quadrupole Response Function of a Trapped Superfluid 
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We derive an explicit relationship between the moment of inertia and the quadrupole response 
function of an interacting gas confined in a harmonic trap. The relationship holds for both Bose 
and Fermi systems and is well suited to reveal the effects of irrotationality of the superfluid motion. 
Recent experimental results on the scissors mode are used to extract the value of the moment of 
inertia of a trapped Bose gas and to point out the deviations from the rigid value due to superfluidity. 
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' Important evidences of superfluidity in Bose-Einstein condensed gases have become available in a recent series of 
experiments. These include the occurrence of quantized vortices the existence of a critical velocity for the motion 
of external impurities |^ , and the rotational properties associated with the scissors mode Q . In particular this latter 
experiment has revealed the consequences of the quenching of the moment of inertia with respect to the classical 
rigid value. The reduction of the moment of inertia follows from irrotationality which makes it impossible, for a 
, superfluid, to respond to a transverse probe and to rotate in a rigid way. The moment of inertia characterizes the 
response of the system to a rotational field —flJz and is defined by the ratio 
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, between the angular momentum (J^) induced by the rotation (hereafter assumed to take place around the z-axis) 
£^ ■ and the angular velocity fl. In superfluid helium the moment of inertia was determined long time ago by rotating the 
bucket containing the liquid and by measuring the angular momentum of the sample ||^. This experiment pointed 
out the occurrence of signiflcant deviations from the rigid value for temperatures below the A-point. Relation (0) 
' naturally defines the moment of inertia in the linear regime. At higher angular velocities new phenomena, associated 
with the creation of quantized vortices, take place in superfluids M. Like in superfluid helium, also in a cold trapped 
gas the moment of inertia is expected to deviate significantly from the rigid value. The temperature dependence of 
^ , O for a trapped Bose gas below the critical temperature Tc for Bose-Einstein condensation and the corresponding 
I ' deviations from the rigid value were investigated in |^ . First estimates of the moment of inertia of a trapped Fermi 
■ O ■ gas below the BCS transition have been also reported In a trapped gas one can rotate the confining potential 
^ using magnetic or optical methods (l^J^ . Actually this procedure has proven quite successful in producing quantized 
^ ^ vortices in Bose-Einstein condensed gases However, the direct measurement of the angular momentum induced 
• • , by the rotation of the trap is difficult in atomic gases, since most of diagnostic techniques, based on optical imaging, 
. ^ ' provide information on the density profiles, either in situ or after the expansion of the gas. 

, In this work we propose a method to measure the moment of inertia of a trapped gas based on the study of the 
^ ■ collective oscillations resulting from the coupling between the rotational and the quadrupole motion. Our analysis, 
I which develops the ideas introduced in is based on the derivation of an explicit relationship between the moment 
of inertia and the quadrupole response function, and will be applied to the recent experiments on the scissors mode of 
a trapped Bose gas |^ . The possibility of revealing the effects of superfiuidity using this method would be particularly 
appealing in the case of Fermi systems, where most of diagnostic techniques are insensitive to the occurrence of the 
BCS transition. 

The starting point of our analysis is the commutation rule {H = 1) 

[H, J,] =. -lixVy - yV,)K=xt(r) , (2) 



where Jz = X^iLi ^iPi ~ ViPf the third component of the angular momentum operator and 

N 2 
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is the Hamiltonian of the system, containing the external confining potential Vcxt (r) . Result ^ is independent of the 
choice of the two-body potential, and follows from the rotational symmetry of the internal Hamiltonian H — Vcxt- In 
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the following we will consider harmonic potentials 



Kxt (r) = ^ Y [ujlxj + ulvf + ulzf) (4) 

with an asymmetry (ui^ '^y) in the rotating plane. In this case the commutation rule (H) takes the simple form 

[H,J,]^-zM{lo1^uoI)Q , (5) 

where 

N 

Q^^Xiyi (6) 

is the relevant quadrupole operator. Eqs. (^)-@ explicitly reveal the coupling between the degrees of freedom asso- 
ciated with the angular momentum and quadrupole variables, induced by the deformation of the confining potential. 

For a sample in thermal equilibrium at temperature T the moment of inertia is given, according to linear response 
theory, by the expression 

e ^ 1 y [e-^- - e-^-l MlMl , (7) 
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where \n) and LOn are eigenstates and eigenfrequencies of the Hamiltonian 2^ is the canonical partition function 
and e~'^"" is the usual Boltzmann factor, with /? = l/ZcsT. 

Because of the algebraic relationship (H) the moment of inertia can be also written in the form 

e = — — (w^-cj ) / dcj ^ (8) 



where 



Xq(^) = 2 Yy - e-^"'"]|(m|Q|n)|^,5(^ - c.„„) , (9) 



n,7n 



is the imaginary part of quadrupole dynamic response function 

and w„m = "-^n ~ '-^m- Here and in the following the frequency integrals are taken between — oo and +oo. It is useful 
to compare the moment of inertia (|^) with the rigid value 

e,ig-MiV(x2 + y2) (11) 

given by classical theory. One can easily show that the average (^Tj) is proportional to the double commutator 
of the Hamiltonian (^) with the quadrupole operator (^. In fact the following sum rule holds P: {[Q,[H,Q]]) = 
N{x'^ + y^) /M. This result is the natural extension of the dipole Thomas- Reiche-Kuhn sum rule ]12| to the quadrupole 
operator. On the other hand the double commutator is directly related to the energy weighted moment of the 
imaginary part (H) of the quadrupole response function, so that one finally obtains the exact relationship N{x'^ + y^) = 
—M Jc?wxq (w)i^/7r, which permits to write the ratio between the moment of inertia (H) and the rigid value ( [ll| ) in 
the useful form 

« -l^l-^lf^P^^ (12, 



Eq. ( p^ ) represents a major result of the present paper. It explicitly shows that the occurrence of superfluidity, 
characterized by a reduction of the value of the moment of inertia with respect to the rigid value, must show up 
also in the structure of the quadrupole excitations of the system, characterized by the function Xq{^)- Actually the 
physical consequence of superfluidity is rather clear. The constraint of irrotationality suppresses the modes with low 
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frequency [u) ~ \u]x — As a consequence in a perfect superfluid the integral (||) approaches a constant value when 
LOx LOy with the consequent quenching of the value of 9, which behaves like (w^ — w^)^. The suppression of the 
low frequency oscillations in the presence of Bose-Einstein condensation has been experimentally pointed out in [|j . 
Conversely, in a normal system the occurrence of low frequency modes is crucial to bring the moment of inertia to 
the rigid value. 

It is worth discussing the effects of superfluidity more explicitly. The equations characterizing the macroscopic 
behaviour of a superfluid at zero temperature have the form of irrotational hydrodynamics 
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where ^\oc{p) is the chemical potential of a uniform system evaluated at density p, and v is the irrotational superfluid 
velocity. These equations hold in particular for a cold dilute Bose gas in the Thomas-Fermi limit |l3). In this case 
Mioc(p) = gp is linear in the density, the coupling constant g = ATra/m being fixed by the s-wave scattering length 
a. In this limit the hydrodynamic equations (n3[]lj) are equivalent to the Gross-Pitaevskii equation for the order 
parameter and can be easily solved to obtain the collective frequencies of a trapped Bose-Einstein condensed gas Q . 
In Ref. ||l^ the same Eqs. (p^p^ have been used to calculate the corrections to the collective frequencies of a trapped 
Bose gas, due to the inclusion of the Lee-Huang- Yang term |^ in the equation of state. Applicability of Eqs. ([131-1^) 
is not however restricted to Bose superfluids. They are in fact expected to hold also in the case of Fermi superfluids. 
Of course in this case the equation of state can be very different due Pauli effects. For example, in a very dilute 
Fermi gas the chemical potential is a purely kinetic effect: piodp) = {Sn'^ p)'^^^ /2m (we consider here a Fermi gas of 
two spin species with equal density). In this limit the hydrodynamic equations of superfluids have been investigated 
in |l^. In general the applicability of Eqs. (p^p^ is subject to the following constraints: (i) the non-superfluid 
component of the system should be negligible and hence the temperature should be much smaller than the transition 
temperature to the normal phase; (ii) the sample should be large enough in order to ensure the validity of the local 
density approximation for the chemical potential (Thomas- Fermi limit); (iii) the length scale of the oscillations should 
be larger than the healing length. 

Starting from the equations of irrotational hydrodynamics it is straightforward to calculate the moment of inertia. 
To this purpose it is convenient to rewrite these equations in the frame rotating with angular velocity il. One finds 



|^ + V(p(v-J^ xr)) =0 
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The stationary solutions in the rotating frame have the simple irrotational form 



V{xy) 



(17) 



From Eq. (17) one can easily calculate the angular momentum and, using definition (|^), one finds the irrotational 
result pi: 
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(18) 



for the moment of inertia, showing that in a superfluid the value of the moment of inertia is smaller than the rigid 
value ||2^]. Results ( p7[ - p^ ) are not restricted to the perturbative regime of small il. They are simply the consequence 
of the fact that the equilibrium density, fixed by the stationary solution of the Euler equation dlq) , depends on the 
variables x and y through the quadratic combination a^x^ -\- Uyy^ [ p^ . For small angular velocities the moment of 
inertia ( p^ ) reduces to 



e = 



H 2 



Orig 



(19) 



Another important feature emerging from the hydrodynamic Eqs. ( |13[jl4D is that also the behaviour of the surface 
oscillations, characterized by the condition V • v = 0, is independent of the equation of state |21|. In particular, by 
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adding a term of the form Xxye^'^* + c.c. in the external potential one finds that the quadrupole operator (g) excites, 
in the linear limit, only one mode with frequency 



t^HD = y ^2 + tj2 ^ (20) 

SO that the imaginary part of the quadrupole response function takes the simple form 



By inserting result (|2l|) into Eq. one immediately recovers the irrotational value (|19|) for the moment of inertia. 

In the second part of the work we show that the deviations of the moment of inertia from the rigid value can be 
measured through a rather straightforward experiment which consists of a sudden rotation of the confining trap by 
a small angle, followed by the imaging of the oscillations of the gas. This experiment has been recently carried out 
Q in a trapped Bose gas where the effects of superfluidity predicted in for the scissors mode have been directly 
observed. 

The abrupt rotation of the harmonic trap is characterized by the transformation of coordinates x ^ x cos (po + 
ysm.(j)Q, y — > —xsmcjiQ + ycos(j)Q. For small values of (/)o, satisfying the condition (/)o <C (w^ — w^)/(a;^ + w^), the 
corresponding change in the harmonic potential (||) produces a time dependent perturbation of the form 

i?pcrt = \Qe{t) , (22) 

where A — M{uj'^ — w^)0o, Q is the quadrupole operator ^ and d{t) is the usual step function. Eq. ( p^ can be 
rewritten in its Fourier components as 

^fpert = Ag— / duj^— + B.C. . (23) 
in J uj + If] 

According to linear response theory the average value Q{t) = {Q) of the quadrupole operator evolves, for t > 0, as 



Q{t) = A— / du——XQiio) + c.c. , (24) 
Att J oj + trj 

where xq{^) is the quadrupole response ([lO|). Integrating ( p^ ) in the complex plane and using the analytic properties 
of xq(w), one can rewrite Eq. (|2^ ) in the form 

Q{t) = J dujQiuj)e''^' - J dwQ{io) , (25) 



with 



OH = . (26) 

TT LO 



Combining the above results with Eq. one finally finds the useful relationship 



Origid ^ I dLjQ{iLj)Lj'^ 

which provides the sought link between the moment of inertia of the system and the Fourier transform of the observable 
quadrupole signal 

0(t)- Jdrpir,t)xy . (28) 

In actual experiments the quantity which is more easily measured is the angle of rotation (p{t) rather than the 
quadrupole moment (^8|). The angle (/)(t) fixes the position of the symmetry axis of the gas with respect to the new 
symmetry axis of the confining trap. During the rotation the density of the gas varies according to the law, valid for 
small angles, 

p(r, t) - po{r) + {xVy - yV,)po(r)((/'(t) - M , (29) 
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where /Oo(r) is the density before the sudden rotation of the trap. Using Eqs. ( |2q ) and ( p9| ) one derives the useful 
relationship 

Q{t)^N{x^ -y^){m~M (30) 

between the angle of rotation and the quadrupole moment. 

We are now ready to use our formalism to extract the moment of inertia from the experiment of Ref. [Q. At 
low temperature the angle <j){t) was measured at different times after the sudden rotation of the confining potential. 
The experimental results have confirmed the occurrence of a single frequency in 4>{t), as predicted by hydrodynamic 
theory: 

(j}{t) = 00 cos(wHDi) ■ (31) 

The observed frequency uo coincides, with excellent accuracy, with the hydrodynamic value (po|). Using the relationship 
( ^ ) and inserting the experimental values ojx = 27r(90±0.2) sec^^, ujy — 27r(249±0.6) sec~^ and w = 27r(265.6 ± 0.8) 
sec^^ into Eq. (P7|), one obtains the value 0/6rigid — 0.6, providing an unambiguous proof of the superfluidity of the 
sample. 

At higher temperature, well above Tc, the situation changes drastically. In first approximation the effects of 
interactions can be neglected and the time evolution of (f>{t) can be described by the non-interacting model 

= ^00 [cos {uu+t) + cos {uj-t)] , (32) 

where 

LU± = \uJj: ± LUy\ (33) 

are the two frequencies excited by the external perturbation. The experimental results of |4|, carried out at T '--^ STcj 
well agree with these predictions. In fact the observed frequencies w+ = 27r(338.5±0.8) sec and uj- = 27r(159.1±0.8) 
sec~^ are very close to the theoretical values (^3|) and enter the measured signal (j){t) with practically equal weight. 
These results, inserted into Eq. (p7|), give finally 6/0rigid — 1, confirming the crucial role played by the mode with 
the low frequency co — \ujx ~ i-^yl in providing the rigid response to the rotational field. 

A similar behaviour is expected to occur also in a Fermi gas. In fact at temperatures much smaller that the BCS 
transition the response is still given by the hydrodynamic predictions (|2^) and (pl|). On the other hand, in the normal 
phase of a very dilute Fermi gas one can use the predictions (p^-(|33|) of the non-interacting model. 

The above procedure can be also used to investigate the effects of superfluidity as a function of temperature. 



In particular, by investigating the measured signal ( pq ) and inserting its Fourier transform into Eq. (27), one can 
determine the value of the moment of inertia at any value of T. One should however take into account that, below 
the critical temperature, the density profile of the oscillating system cannot be simply parametrized through a single 
angle of rotation, like in Eq. (p9|), because the normal and the superfluid components will respond in a different way 
to the rotation of the trap. As a consequence one should introduce two distinct angles and the density of the gas will 
oscillate, in the linear limit, according to 

p(r, t) = po{r) + {xVy - y\7x)ps{r)(j)s{t) + {xVy - yV x)pN{r)4>N{t) , (34) 

where po(r) = Ps{t^) + Pn{y) is the sum of the superfluid and normal component at equilibrium. The careful 
investigation of the time evolution of the density profile would provide unique information on the shape of the 
superfluid and normal components. In a dilute Bose gas these components are expected to coincide, respectively, with 
the condensate and the thermal cloud, except for temperature smaller than the chemical potential. The situation is 
very different in Fermi gases, where the total density profile is scarcely affected by the BCS transition |2^, and the 
separation between the two components is much less trivial. 

In conclusion we have shown that the moment of inertia O of a gas trapped by harmonic potential can be explicitly 
measured by analyzing the quadrupole oscillations generated by the a sudden rotation of the trap. An explicit 
evaluation of 8 has been provided using the recent data on the scissors mode available for a trapped Bose gas either 
at very low temperature and above the critical temperature. In the first case the method explicitly reveals the 
quenching of the moment of inertia due to superfluidity. The same method can be naturally applied to a Fermi 
system in order to reveal the effects of superfluidity associated with the BCS transition. Moreover, with suitable 
modifications, it can be also extended to the case of non-harmonic trapping. 

Stimulating discussions with D. Guery-Odelin are acknowledged. We are also grateful to O. Marago for providing 
us the data of the Oxford experiment on the scissor mode. 
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